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A relativistic distorted-wave impulse-approximation model to quasielastic neutral-current 
neutrino-nucleus scattering is developed, using the relativistic mean field theory for the 
bound state and taking into account the final state interaction in the relativistic scattering 
state. Results for the neutrino (antineutrino) reaction on 12 C target nucleus are presented 
| in an energy range between 500 and 1000 MeV. The sensitivity to the strange quark content 

f^S ■ of the axial-vector form factor is discussed. 
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INTRODUCTION 



The neutral-current scattering of neutrinos and antineutrinos on nucleons and nuclei has gained 
in recent years a wide interest in order to determine the structure of hadronic weak neutral current. 
In a measurement of the proton spin structure function, the European Muon Collaboration 0] found 
out a disagreement with the Ellis- Jaffe sum rule Q], when only up and down quark and antiquark 
contribution to the proton spin was taken into account, thus indicating that other flavors might 
contribute to the spin structure of the proton. A measurement of neutrino (antineutrino)-proton 
elastic scattering at Brookhaven National Laboratory (BNL) 0] supported this result by reporting a 
non-zero value for the strange axial- vector form factor of the nucleon. However, the BNL data were 
, critically reanalyzed in Ref. 0], where also strange vector form factors were taken into account. 

It turned out that the BNL data cannot provide us decisive conclusions about the strange form 
factors 0| . The BooNE experiment [j| at Fermi National Laboratory (FermiLab) aims, through 
the FermiLab Intense Neutrino Scattering Scintillator Experiment (FINeSSE) [jj], at performing a 
detailed investigation of the strangeness contribution to the proton spin via neutral current elastic 
scattering. 

^ In order to reach a definite conclusion, nuclear structure effects have to be clearly understood. 

Quasielastic electron scattering calculations 0], which were able to successfully describe a wide 
number of experimental data, can provide us a useful tool to study neutrino-induced processes. 
However, a careful analysis of neutrino- nucleus reactions introduces additional complications, as 
the final neutrino cannot in practice be measured and a final hadronic probe has to be detected: 
the corresponding cross sections are therefore semi-inclusive in the hadronic sector and inclusive 
in the leptonic one. 

General review papers about neutrino-nucleus interactions can be found in Refs. [1,0,0,0]. 
The relativistic Fermi gas (RFG) model was applied to study neutral-current neutrino-induced 
knockout from nuclei in Ref. |l2j], where binding energies corrections and strange-quark axial- 
vector and vector effects were investigated. Detailed analyses of nuclear structure effects on the 
determination of strangeness contribution were performed in Refs. [lil. Il4|, In addition to the 
RFG model a "hybrid" model in which bound nucleons are described by harmonic oscillator wave- 



functions was proposed in Ref. 13], while a relativistic shell model was used in Ref. [14]. The 
effects of final state interactions (FSI) were also studied in Ref. [3] in the framework of random 
phase approximation (RPA) and a continuum RPA model was developed in Ref. ^(|, where 
nucleosynthesis processes were also discussed. The FSI effects are generally large on the cross 
section, but they are reduced when studying the ratio of proton to neutron cross sections. This 
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was discussed in Refs. [l2L Il7l| . The effects of two-body relativistic meson exchange currents 
in neutrino-nucleus scattering at low and intermediate energies were evaluated in Refs. [isL Hfll]. 
The sensitivity of the neutrino-nucleus cross sections to the strange-quark contribution was also 
examined in a relativistic plane wave impulse approximation (RPWIA) in Ref. jijj, where a 
model-independent approach based on eight structure functions was developed. 

In this paper we present a relativistic distorted-wave impulse-approximation (RDWIA) calcu- 
lation of neutral-current v- and P-nucleus reactions in the quasielastic region, where a neutrino 
interacts with only one nucleon in the nucleus and the others nucleons behave as spectators. The 



RDWIA treatment is the same as in Refs. |21l . VIA . |23|, where it was applied to electromagnetic 
knockout reactions and in Ref. where a relativistic Green's function approach was developed 
dealing with the quasielastic charged-current i/-nucleus scattering. The relativistic bound state 
wave-functions are solutions of a Dirac equation contai ning scalar and vector potentials obtained 
in the framework of the relativistic mean field theory \2nL l26j |. The FSI are taken into account 
through a relativistic optical model potential 



27, urn. The effective Pauli reduction has been 



adopted in this case and the resulting Schrodinger-like equations are solved for each partial wave. 

The formalism is outlined in Sec. HU Results are presented and discussed in Sec. IIII1 Some 
conclusions are drawn in Sec. IIVI 



II. THE FORMALISM OF THE SEMI-INCLUSIVE QUASIELASTIC SCATTERING 

The neutral-current ^(i^)-nucleus cross section for the semi-inclusive process can be written as 
the contraction between the lepton tensor and the hadron tensor, i.e., 

da = -2nL» ^pp, (1) 

where G ~ 1.16639 x 10~ n MeV~ 2 is the Fermi constant, and kf = (e^fej), k^ = (e, k) are 
the four-momentum of the incident and final neutrino or antineutrino, respectively, and p N is the 
momentum of the emitted nucleon. 

The lepton tensor has a similar structure as in electromagnetic knockout or in charged-current 
z/-nucleus scattering and can be written as in Ref. [t1. Eol. E3| . After projecting into the initial and 
the final neutrino (antineutrino) states, it is decomposed into a symmetrical and an antisymmetrical 
component, i.e., 

IT = — [igr T IT] , (2) 



E4.E 



with 



I £i — k^ k — f— k^ k^ k% ' k 

I A = i ^ ^ki a kp, (3) 

where e^ ua @ is the antisymmetric tensor with €0123 = — e 0123 = 1. The upper (lower) sign in Eq. [21 
refers to v{p) scattering. 

Here, we assume the reference frame where the z-axis is parallel to the momentum transfer 
q = ki — k and the y-axis is parallel to k{ x k. 

The hadron tensor for the semi-inclusive scattering is given in its general form by bilinear 
combinations of the transition matrix elements of the nuclear weak neutral-current operator 
between the initial state |^o) of the nucleus, of energy Eq, and the final states of energy Ef, given 
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by the product of a discrete (or continuum) state |n) of the residual nucleus and a scattering state 
Xp of the emitted nucleon, with momentum p N and energy E^. It can be written as 

q) = J>; X g I J"(g) | * )(*0 I J 14 (g) I n; X g) 5(E + u - Ef) , (4) 

n 

where the sum runs over all the states of the residual nucleus. 

Here, the transition matrix elements are calculated in the first order perturbation theory and 
in the impulse approximation, i.e., the incident neutrino interacts with only one nucleon, while 
the other ones behave as spectators. Thus, the transition amplitude is assumed to be adequately 
described as the sum of terms similar to those appearing in the electron scattering case |2l| 



in; x g | J"(g) | * > = {xfe? I f(q) I <Pn) , (5) 



where cp n = (n|\Po) describes the overlap between the initial nuclear state and the final state of the 
residual nucleus, corresponding to one hole in the ground state of the target. The single-particle 
current operator related to the weak neutral current is 

f = F^Q 2 )^ + *^2 V (Q W - Ga(QV7 5 + iMQ V7 5 , (6) 

where k is the anomalous part of the magnetic moment, = (ui,q), Q 2 = \q\ 2 — to 2 , is the 
four-momentum transfer, and a^ v = (i/2) [7^,7"]. F^ and are the isovector Dirac and Pauli 
nucleon form factors, taken from Ref. |29|. Ga is the axial form factor and Fp is the induced 
pseudoscalar form factor, which does not contribute to massless neutrino scattering. The vector 
form factors F^ can be expressed in terms of the corresponding electromagnetic form factors for 
protons (Ff) and neutrons (Ff), plus a possible isoscalar strange-quark contribution (Ff), i.e., 

F7 = Q - 2 sin 2 # w ) Ff - l -F? - ~F? (proton knockout) , 
pV _ f 1 n^2o\ T?n 1 kip 1 r s 
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2 sin 2 6» w ^ Ff - ^Ff - ]^Ff (neutron knockout) , (7) 



where #w is the Weinberg angle ( sin 2 6w — 0.2313). Standard dipole parametrization is adopted 
for the vector form factors [l^, L^l • The axial form factor is expressed as [2!j] 

G A = ^ (9 A ~ 9a) g (proton knockout) , 
1 

Ga = — - (pa + 9%) G (neutron knockout) , (8) 

where g\ ~ 1.26, g s A describes possible strange-quark contributions, and G = (1 + Q 2 /M]f)~ 2 with 
Ma = 1.026 GeV. The strange vector form factors are taken as [15( 

- F 2 S (0) 



F!(Q Z ) 



(l + r)(l + QVM 2 )2' 



FlS(Q2) " (l + r)(l+ Q QVM 2 )2' (9) 

where r = Q 2 /4M 2 , iff (0) = Ms, = ~(r 2 )/Q, and M v = 0.843 GeV. The quantity u s is the 
strange magnetic moment and (r 2 ) the squared "strange radius" of the nucleus. 
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The differential cross section for the quasielastic neutral-current ^(^)-nucleus scattering is ob- 
tained from the contraction between the lepton and hadron tensors, as in Ref. Q. After performing 
an integration over the solid angle of the final nucleon, we have 



da 



dedQ dT\ 



G 2 2 2 $ 
2^ £ C ° S 2 



V R 00 + V ZZ R ZZ - V 0z R 0z + V T Rt ± ^y-Rxj/ 



(2vr) 3 



(10) 



where $ is the lepton scattering angle and the relativistic energy of the outgoing nucleon. The 
coefficients v are given as 



vo 



1 



voz = ,— r 



2 ^ Q 2 
UT = tan 2 + 2|^f 



tan ■ 



, 2® ,Q 2 

tan — + - — - 



(11) 



where the neutrino mass has been neglected. 

The response functions R are given in terms of the components of the hadron tensor as 



Rzz 

Roz 
R T 



J d^ N W 00 , 
dft N W zz , 
dS^ N 2 Re{W 0z ) , 
dft N (W xx + W yy ) 
dft N 2 Im(W^) . 



(12) 



Here, we are interested in the single differential cross section with respect to the outgoing 
nucleon kinetic energy Tn, i.e. 



da 
dTV 



da 



de dn dT N 



dedO 



and in the "total" cross section 



a 



da 



dT N . 



(13) 



(14) 



In the calculation of the transition amplitudes of Eq. Elthe relativistic single-particle scattering 
wave function is written as in Refs. 2l|, |30j in terms of its upper component, following the direct 
Pauli reduction scheme, i.e., 



X Pn 



f+ 



M + E + S\E) - V\E) 



(15) 



where S{E) and V(E) are the scalar and vector energy-dependent components of the relativistic 
optical potential for a nucleon with energy E ji^]. The upper component, \&f + , is related to a 
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two-component spinor, <3?f, which solves a Schrodinger-like equation containing equivalent central 
and spin-orbit potentials and which is obtained from the relativistic scalar and vector potentials 



28, 3 



i.e. 



f+ 



D(E) 



1 + 



D1(E) $ f , 
S(E) - V{E) 



M + E 



(16) 



where D(E) is the Darwin factor. 

The single-particle overlap functions (p n in Eq. |S] are taken as the Dirac-Hartree solutions of a 
relativistic Lagrangian containing scalar and vector potentials ijj 2f|. 



III. THE RESULTS FOR 12 C 

The calculations have been performed for the 12 C nucleus with the same bound state wave 
functions and optical potentials as in Refs. [ill E3. liil. liil. liol] . where the relativistic approach was 
developed to study (z^,/x~), (e,e'p), (7,p), and (e,e') reactions. 

The relativistic bound state wave functions have been obtained from Ref. jij|, where rela- 
tivistic Hartree-Bogoliubov equations are solved in the context of a relativistic mean-field theory 



and reproduce single-particle properties of several spherical and deformed nuclei 26j. The scat- 



tering state is computed by means of the energy-dependent and A-dependent EDAD1 complex 



phenomenological optical potential of Ref. |27( , which is fitted to proton elastic scattering data on 
several nuclei in an energy range up to 1040 MeV. The initial states (p n are taken as single-particle 
one-hole states in the target, with a unitary spectral strength. The sum runs over all the occupied 
states. In this way we include the contributions of all the nucleons in the nucleus, but disregard 
effects of correlations. These effects, however, are expected to be small on the semi-inclusive cross 
section, and, moreover, should not change the effects of FSI and of the strange-quark content of 
the form factors, which represent the main aim of the present investigation. 

Calculations have been performed for three different values of the incoming energy, i.e., 
F v (v) = 500, 700, and 1000 MeV. Firstly, we investigate the effects of the FSI between the emitted 
nucleon and the residual nucleus, which are described in our approach by means of a phenomenolog- 
ical optical-model potential. The cross sections calculated in RDWIA and in RPWIA are compared 
in Fig. ^for proton emission and in Fig. [2] for neutron emission. A reduction of the cross sections 
of ~ 50% is found in Fig. ^ similar to the one obtained in the case of electromagnetic proton 
knockout 7]. However, a larger reduction (~ 60%) is obtained in Fig. |2]when a neutron is emitted. 
Therefore, the ratio between the cross sections for proton and neutron emission is enhanced by 
distortion. The comparison between the cross sections for an incident neutrino and antineutrino 
gives the contribution of the axial component of the neutral current. The results are in reasonable 
agreement with those of Ref. jjfl ] . 

The effect of a strange-quark contribution to the axial-vector form factor is studied in Figs. 
El and 0J where the cross sections with g\ = —0.19 are compared with the corresponding ones 
without any strange-quark content in the weak current. The results with g^ = —0.19 are enhanced 
by ~ 50% with respect to those with = in the case of proton knockout and reduced by 
~ 30% in the case of neutron knockout. This different behavior is due to the different sign of the 
interference term gA9% in the proton and in the neutron case. The effect of g^ is less significant in 
the cross sections with an incident antineutrino. 

The effect of the strange vector form factors F^Q 2 ) and Fl(Q 2 ) on the neutrino quasielastic 
cross sections is shown in Fig. [SJ We have chosen as in Ref. jj] = —0.40 and Ff = — (r 2 )/Q = 
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0.53 GeV -2 . The Q 2 dependence is given in Eq. |UJ While FJ decreases the cross sections, Ff gives 
an enhancement that, at E v = 1 GeV, almost cancels the effect of Ff . These effects are smaller 
than that produced by the strange axial-vector form factor, g^, which is shown also in Fig. [21 for 
a comparison. 

In order to separate the effects of the strange-quark contribution and of FSI on the cross sections, 
it was suggested in Refs. 2) 13- HI to measure the ratio of proton to neutron yields, as this ratio is 



expected to be less sensitive to distortion effects than the cross sections themselves. Moreover, from 
the experimental point of view the ratio is less sensitive to the uncertainties in the determination 
of the incident neutrino flux. In Fig. EJthe ratio for an incident neutrino is displayed as a function 
of the outgoing nucleon kinetic energy. It is very sensitive to and exhibits, for g^ = —0.19, 
a pronounced maximum at T/v — 0.6 E v . The RPWIA results are similar to those of Ref. |2oT |. 
A sizable enhancement of the ratio is produced by FSI. This result is mainly due to the different 
effects of the optical potential in proton and neutron knockout. It means that the argument of 
looking for the strange-quark content in this ratio is strengthened by distortion. 

The global effect of the strange-quark contribution is shown in Fig. where the cross sections, 
integrated over the emitted nucleon kinetic energy, are displayed. The cross section is increased by 
~ 50% for proton knockout and decreased by ~ 30% for neutron knockout. The effect is smaller 
for the cross sections with an incident antineutrino. The ratio of proton to neutron cross sections is 
shown in Fig. El It is largely enhanced (by a factor ~ 2) for neutrino scattering, and less enhanced 
(~ 50%) for antineutrino scattering. 

Finally, we compare our results with the data of the BNL 734 experiment Wide band 
neutrino and antineutrino beams of average energies of 1.3 and 1.2 GeV, respectively, were used 
to study neutral-current z/(z/)-proton scattering. Approximately 80% of the events were due to 
quasielastic scattering from 12 C nuclei and 20% to elastic scattering from free protons. Experimen- 
tal results were presented in the form of a flux-averaged differential cross section per momentum 
transfer squared Q 2 , which, for a free proton, is related to the recoil energy by the simple relation 
Q 2 = 2MT p . In the case of quasielastic scattering on a complex nucleus, this is only approximately 
true. However, we can still define 

O = Wl I ^{T„ m )^r , (17) 

where the v(p) spectrum, $(T v r v -\), is normalized to 1. In Fig. |5]our results for the flux-averaged 
cross section are displayed and compared with BNL data. The quasielastic cross sections have been 
weighed with the free proton scattering ones 0, E3l- We show results with g s A = —0.19 and without 
the strange-quark contribution. Moreover, following Ref. 0|, we give the effect of including the 
strange vector form factors, Ff and Ff, with Ff = — (r 2 )/6 = 0.53 GeV -2 , Ff(0) = —0.40, and the 
Q 2 dependence given in Eq. |HJ The strange-quark contribution produces an enhancement of the 
cross sections, which makes them slightly higher than the experimental data. The strange weak 
magnetic contribution decreases the cross section, while the axial and weak electric components give 
an enhancement. In the lower panel of Fig. El the ratio of neutrino to antineutrino flux-averaged 
cross sections is presented and compared with the BNL data. The shape of the experimental data 
is reasonably reproduced. 



IV. SUMMARY AND CONCLUSIONS 

We have presented relativistic calculations for neutral current z/(z^) -nucleus quasielastic scatter- 
ing. The reaction mechanism is assumed to be a direct one: the incident neutrino (antineutrino) 
interacts with only one nucleon in the target nucleus and the other nucleons remain as spectators. 
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A sum over all single particle occupied states is performed, using an independent particle model 
to describe the structure of the nucleus. The scattering state is an optical-model wave function. 
Results for the 12 C target nucleus have been presented at neutrino (antineutrino) energies up to 
1000 MeV. The optical potential produces a large reduction of the cross sections. The sensitivity to 
the strange-quark content of the axial-vector form factor has been investigated. A negative value 
of = —0.19 sensibly increases the cross section for proton knockout and decreases it for neutron 
knockout. A large enhancement of the ratio of proton-to-neutron cross sections is obtained with 
respect to the case with <^ = 0. The enhancement is increased by FSI. The sensitivity to the 
strange-quark contribution of the vector form factors has also been discussed. FJ decreases and 
Ff enhances the cross sections. The two effects tend to cancel each other and are anyhow smaller 
than the effects produced by a strange axial-vector form factor. Finally, a comparison with the 
BNL 734 experimental data has been presented. The shape of the experimental data is reasonably 
reproduced. 
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FIG. 1: Differential cross sections of the v{p) quasielastic scattering on 12 C as a function of the outgoing 
proton kinetic energy T p . Solid and dashed lines are the results in RDWIA and RPWIA, respectively, for 
an incident neutrino. Dot-dashed and dotted lines are the results in RDWIA and RPWIA, respectively, for 
an incident antineutrino. The strangeness contribution is here neglected. 
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FIG. 2: The same as in Fig. ^but for neutron knockout. 
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FIG. 3: Differential cross section in RDWIA of the v{p) quasielastic scattering on 12 C as a function of the 
outgoing proton kinetic energy T p . Solid and dashed lines are the results with g\ = and g\ = —0.19 in 
the case of an incident neutrino. Dot-dashed and dotted lines are the corresponding results for an incident 
antineutrino. 



10 




100 200 J00 400 500 



T n [MeV] 

FIG. 4: The same as in Fig. OUbut for neutron knockout. 
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FIG. 5: Differential cross section of the v quasielastic scattering on 12 C as a function of the outgoing proton 
kinetic energy T p . Solid lines are the results with no strangeness contribution, dashed lines with g\ = —0.19, 
dot-dashed lines with g% = -0.19 and Ff (0) = -0.40, dotted lines with g% = -0.19, Ff (0) = -0.40 and 
Fl = -(r s 2 )/6 = 0.53 GeV" 2 . 
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FIG. 6: Ratio of proton to neutron cross sections of the v quasielastic scattering on 12 C as a function of 
the outgoing nucleon kinetic energy Tn- Solid and dashed lines are the results in RDWIA with — and 
g\ = —0.19. Dot-dashed and dotted lines are the same results but in PWIA. 




E v [MeV] 

FIG. 7: The total cross section of the v(p) quasielastic scattering on 12 C integrated over the outgoing nucleon 
kinetic energy. Upper panel: proton knockout. Lower panel: neutron knockout. Line convention as in Fig. 

m 



12 




tMeV] 



FIG. 8: Ratio of proton to neutron total cross sections of the v[y) quasielastic scattering on C as a function 
of the incident neutrino (antineutrino) energy. Solid and dashed lines are the results in RDWIA with g s A = 
and g\ = —0.19. Dot-dashed and dotted lines are the same results but for an incident antineutrino. 
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FIG. 9: Upper panel: differential cross sections of the v{v) quasielastic scattering, flux-averaged over BNL 
spectrum as a function of the momentum transfer squared. The four upper curves are for incident 
neutrino and the four lower ones for incident antineutrino. Lower panel: ratio of neutrino to antineutrino 
flux-averaged cross sections. Line convention for both panels as in Fig. |S| Experimental data from Ref. Q. 



